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ABSTRACT 

Context. The impact of quasars on their galaxy neighbours is an important factor in the understanding of galaxy evolution models. 
Aims. The aim of this work is to characterize the intermediate-scale environments of quasars at low redshift (z < 0.2) with the most 
statistically complete sample to date using the seventh data release of the Sloan Digital Sky Survey. 

Methods. We have used 305 quasar-galaxy associations with spectroscopically measured redshifts within the projected distance range 
of 350 kpc, to calculate how surface densities of galaxies, colors, degree of ionization, dust extinction and star-formation rates change 
as a function of the distance to our quasar sample. We also identify the companion Active Galactic Nuclei from our main galaxy 
sample and calculate surface density for dilferent galaxy types. We have done this in three different quasar-galaxy redshift difference 
ranges |Az| < 0.001, 0.006, and 0.012. 

Results. Our results suggest that there is a significant increase of the surface density of blue neighbours around our low-redshift 
quasar sample that is steeper than around non-active field galaxies of the same luminosity and redshift range. This may indicate that 
quasar formation is accomplished via a merging scenario. No significant changes in star formation rate, dust extinction, degree of 
ionization or color as a function of distance from the quasars was observed. We could not observe any direct effects from quasars on 
the their companion galaxies. 

Key words, active - nuclei - star formation - quasars: general - interactions - evolution 



1. Introduction 



Since their discovery in 1963 (ISchmidtl 1 19631) quasars have 
sparked widescale interest among astronomers. Their large lu- 
minosities and variation over redshift are incomparable to that of 
any other galaxy type, while their decline in number in our Local 
, Universe remains poorly understood. Low-redshift quasars are 
' not only under-represented, but also considerably fainter com- 
' pared to their high-redshift members. Regardless of whether they 
were victims of their own high-energy consuming nature, or 
■ rather strangled their less luminous surrounding galaxies, traces 
I of their dynamic history could be found by observing in which 
i way these most extreme active galactic nuclei influence the evo- 
' lution of their galaxy neighbours. These effects could mark ei- 
ther direct influence of AGN, or trace the environments most 
favorable to their occurrence. 

The extremely luminous nature of quasars is believed to be 
driven by an active galactic nucleus due to accretion of mate- 
rial upon a supermassive black hole (e.g. iLvnden-Belll 119691: 
lReesLfT98l . Feeding such a black ho le may affect the formation 
of ga laxies in proximity to quasars (iGunni 119791: ISilk & ReesL 
Il998 l). and this could play a large role in structure formation 
scenarios derived from standard cosmological models. 

Could quasars be the result of mergers according to a hi- 
erarchical structure formation scenario predicted by Cold Dark 
Matter Models? For many years there have been indications 
that both nuclear activity and star formation may be triggered 
or enhanced by interactions between galaxies. For example. 



the presence o f post-sta rburst popula tions i n close quasar com- 
panio ns (e.g. Canalizo & Stocktonl 119971: [Scannapieco et all 
2005t lAlmeida et alT 1201 Ih . morphologic al asymmetries with 
high star formation rates near quasar hosts (iGehren et al.l 119841 : 
Hutchings et al., 1984) and possible quasar companion galax- 
ies (Can alizo & Stocktonl 1200 1|). Ion i zed g as has also been ob- 
served around quasars (e.g. iBorosonI 119851) . while others stud- 
ies have suggested that high-luminosity quasars cou ld ionize the 
galactic medium up to several megaparsecs away (iReesl 119881: 
Babul & White, 1991). 

A number of additional studies have suggested that quasars 
influence the evolution of normal galaxies. There are essentially 
two different populations of normal "non-active" galaxies; red 
mostly quiescent galax ies with early- t ype m orphology, a nd blue 
star- f orming galaxies (IStrateva et al.l l2001t: .Brinchmann et all 
l2004l) . The color bimodality of these galaxy populations sug- 
gest that they originate from a mechanism that quenches star 
formation in the blue galaxy population turning them red. This 
hypothesis has resulted in a number of ideas of possible feed- 
back mechanism s associate d with a ctive gala ctic nuclei (e.g . 
Tabor &Binrievl [1993: HopkinseLalJ, 2006; Sii acki & Springej 



2006"). These feedback mechanisms have been proposed as an 
alternative to more conservative ex planations such as merg- 
ers between late-type gala xies (e.g. iToomre & Toomrel 1 19721 
Barnes & Hemauist', '1992") or loss of gas due to ram pressure 
stripping (Gunn & Gott, 1972). 

However, more study is needed in order to distinguish the 
various scenarios presented above. A great deal of work has 
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gone int o the study of close pairs (less than 30 kpc) of galax- 
ies (e.g. iRogers et al.L l2009 l) and of the surrounding environ- 
ments of quasars at larger radii (larger than 1 Mpc) most re- 
cendy by Coldwell & Lambas (2006) hereafter CW06, but there 
have been few detailed studies on an intermediate scale. In the 
previous larger-scale study of quasar environments up to 3 Mpc 
by CW06, a small (20%) excess of blue star-forming galaxies 
around quasars was shown for redshifts z < 0.2. CW06 also 
demonstrated that quasars tend to avoid high-density environ- 
ments such as clusters and galaxy grou ps. This conclu sion has 
been supported by other studies (e.g. iLietzeni |2009|) . Active 
Galactic Nuclei ( AGN) also tend to avoid high-d ensity environ- 
ments (e.g. Way e t~anil998l:lAmold et al.Ll2009h . 

Many studies of quasar environments using the Sloan Digital 
Sky Survey (SDSS) have used galaxies with photometrically de- 
termined redshifts (photometric redshifts, Photo-Z) which have 
large redshift uncertainties 6~ > 0.025. These allowed more 
quantitative studies than qualitative, which underlines the im- 
portance of careful quasar-galaxy association studies with spec- 
troscopically defined redshifts with small errors 5, < 0.001. 

In the CW06 investigation SDSS quasars were restricted to 
redshifts z < 0.2 and luminosities brighter than M^ < -23. The 
aim was to explore the distributions of quasars in different envi- 
ronments. They also examined the distribution of the quasar col- 
ors and magnitudes, as well as the morphology and spectral type 
of the galaxies surrounding the quasars. However, their results 
were difficult to interpret due to the fact that they used fraction of 
total number of galaxies in each histogram bin without inform- 
ing the reader of the total number of galaxies per bin. They also 
relied on the SDSS spectral classification coefficient (eClass) 
to estimate the morphology of the quasar companion galaxies. 
The eClass coefficient is calculated from Principal Component 
Analysis (PCA). We show in Section [TTI that this method is not 
an accurate classification of morphology. 

We use a non-volume-limited sample drawn from the sev- 
enth data release (DR7) of the SDSS. The DR7 contains more 
than 100,000 quasars of which we use a small subsample to 
study how quasars influence their nearby galaxies at low red- 
shifts. We focus on examining how colors, surface densities, star 
formation rates and ionization of the galactic medium of galax- 
ies near quasars change depending on the distance between the 
quasar and the neighbour. This is done for a sample consist- 
ing of 305 quasar-galaxy associations, at redshifts z < 0.2 us- 
ing Standard Cosmology (Qa=0.70, Qm=0.30, H„=70 km s^' 
Mpc"^). In a future publication we plan to replicate this study 
for higher-redshift objects to investigate possible redshift evolu- 
tion in quasar-galaxies associations. With these two studies, it 
may be possible to understand how quasars impact surrounding 
galaxies as a function of time. Such a comprehensive study could 
better constrain AGN feedback models and contribute to the un- 
derstanding of the formation of the most luminous galaxies in 
our Universe. 

In Section|2]we describe our sample selection, Section|3]will 
discuss how we estimate the environment around the quasars. 
Section ID is a discussion of the results of our analysis and con- 
clusions. 

2. Sample Selection 

2.1. Spectroscopic redshift sample 

The quasar and galaxy samples were drawn fr om the SDSS 
dYork et all [2000) DR7 dAbazaiian et al.L l2009l) . The spectro- 
scopic data in this catalogue contains 929,555 galaxies and 



121,363 quasars in five optical bands (u, g, r, i, z) that cover 
the complete wavelength range of the CCD-camera. 

The quasars were extracted from the SDSS quasar catalogue 
(Sc hneider et al.L 12003). These have been selected such that they 
have at least one emission line with full width at half max- 
imum (FWHM) lai-ger than 1000 km s"' and M,- < -22, and 
have highly reliable redshifts (zConf > 0.95). Both the quasars 
and galaxies were chosen to be in the redshift range 0.03 < z 

< 0.2 which means that most of the quasars we have are low 
luminosity quasars and Seyferts rather than conventional high- 
luminosity quasars. The limited apparent luminosity of the spec- 
troscopic selection of quasars implied by Schneider et al. (2002) 
will likely help us to avoid the low luminosity Seyfert Is. Objects 
were excluded from our s ample if they were flagged via the 
SDSS photometric pipeline (Lupton et al.','2001') as overly bright 
(flags&0x2=0), saturatedj;flags&0x40000=()}, or from blended 
images (flags&0x8=0) (iStoughton et al.Ll2002i see Table 9). We 
do a SDSS query for quasars within the redshift range to see 
how many we get with and without the saturation, blending and 
brightness flags. We observed a loss of 5% of the quasars, an ef- 
fect that we consider being negligible. As quasars are the most 
luminous among AGN, we do not expect this to influence the 
number of AGN neighbours later in our analysis. We used the 
BestObjIDs in our queries for both galaxies and quasars, since 
the imaging data in BESTDR7 is of the highest qualitjQ. This is 
particularly important for the careful selection of blue galaxies 
in this study. 

To obtain the quasar neighbours the casjob^ query was done 
in two parts: a low redshift part 0.03 < z < 0. 1 where we searched 
for neighbouring galaxies up to 10' and a high-redshift part 0.1 

< z < 0.2 with companions within 7' of the quasar. We limited 
the redshift difference between the galaxy and the companion to 
be |Az I < 0.012 in the queries so that this may be compared to 
similar Photo-Z studies which have photometric redshift errors 
of this order. The goal is to obtain any galaxy within 350 kpc 
of a quasar in our sample. With the minimum 55" fiber separa- 
tion, pairs closer than 100 kpc from each other at redshifts z <: 
0.1 will not be measured on the same spectroscopic plug plate. 
Many of the pairs with r < 17.7 will be seen in the SDSS imag- 
ing survey, but will lack a redshift measurement for one of the 
objects. For both quasar-galaxy pairs within 100 kpc and with 
redshifts z ^ 0. 1 to be spectroscopically measured they have to 
reside in plate overlap regions. For this reason we have a partic- 
ularly small sample size. 

The following quantities were retrieved for each quasar and 
galaxy in our sample: 

- spectral line (absorption or emission) information 

- apparent dereddened magnitudes (corrected for extinction) 

- K-corrections and rest-frame absolute magnitudes 

- isophotal axis length 

- radii including 50% and 90% of the Petrosian flux in the i- 
band 

We calculated the projected distances for all 305 quasar- 
galaxy associations (there were 284 unique quasars). Galaxies 
with internal extinction corrected colors of Ug- Re < 2.2 (using 
rest-fram e magnitudes) will h enceforth be referred to as "blue" 
galaxies dStrateva et al.Ll2001i) . while the remaining galaxies will 
be called red (see Section [J!2T l. 



' See the SDSS casjobs Schema browser for more information 
http://casjobs.sdss.org/dr7/ en/help/browser/browser.aspl 
" http://casjobs.sdss.org 
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In subsequent calculations of surface densities, star- 
formation rates, ionization degree and colors we have separated 
AGN from normal galaxies (see Section l23t . 

2.2. Line flux and extinction corrections 

The Balmer emission line fluxes and equivalent widths (EW) 
were corrected for underlying stellar absorption by assuming av- 
erage absorption line strengths corresponding to 2.5 A in EW 
for Ha and 4 A for li/3. Internal extinction corrections were only 
done if EW(H/3) > 5 A since Hfi lines lower than 5 A are too 
noisy to permit a reliable extinction correction. The spectral lines 
were internal extinction-corrected following a standard interstel- 
lar extinction curve (e.g. Whitford, 1958; Osterbrock, 1989*). 

The colors were corrected for extinction by taking the incli- 
nation of the observed galaxies into accoun t. The was done for 
all galaxies fulfilling the criteria described in ICho & Parkl(l2009h 
using their derived analytical expressions. Colors were corrected 
for galaxies having < u-r < 4, Ha 0-200 A absolute magnitude 
-21.95 < Mr <-19.95 and a concentration index (CI) in the range 
1.74 <C/< 3.06. 

At earlier stages of our analysis, we did the same study with- 
out any line flux corrections. The addition of line flux corrections 
has not influenced the conclusions. 

2.3. AGN activity 

In some theories AGN can be responsible for inducing star- 
formation activity in nearby galaxies, while in other theories 
they are hypothethized to quench star formation via AGN feed- 
back. In order to investigate either of these possibilities it is nec- 
essary to separate AGN from normal star-forming galaxies in 
the quasar companion sample. In order to separate AGN from 
starforming H ii type galaxies in our neighbour galaxy sample 
we used BPT line-ratios (Baldwin-Phillips & Terl evich. .1981.) 
combined with lKauflinann et al.l (120031) criteria: 

log([Oui]/H/3) > 0. 61 /(logi[N 11] /Ha)) - 0.05) + 1.3 (1) 

To remove the remaining broad-line (type 1 ) AGN an addi- 
tional criterion was added, a-(Ha) > 15 A. With these two crite- 
ria 69 AGN were identified in the companion galaxy sample. 

2.4. Field gaiaxies 

A companion field galaxy sample was obtained from the papers 
of Way et al. (2011, in prep) 16526 galaxies were classified as 
"field" (Iblock=12-19) by the "Bayesian Blocks" method out- 
lined in the aforementioned papers. 16509 of the original "field" 
sample had SDSS casjobs Object IDs. The object IDs were nec- 
essary to obtain additional photometric and spectroscopic pa- 
rameters not present in the original catalogs. Among these we 
selected non-AGN field galaxies using the criteria discussed in 
Section l23] Using the same procedure outlined in Section lCTI to 
obtain neighbours for the quasar sample, we obtained 435 field 
galaxy-galaxy pairs in the redshift range 0.03 < z < 0.1. 

2.5. Photometric redshift samples 

Approximately 1/3 of the Sky is covered by more than one spec- 
troscopic plate, and it is estimated that roughly 2/3 of all pairs 
with angular separation less than 55" will not be detected. This 
may bias our sample towards more wide pairs. To examine how 



the conclusions on eventual clustering around quasars can be in- 
fluenced by the fiber collision constraint, we constructed a pho- 
tometric redshift sample with neighbours from the Photo-z cat- 
alogue. For the quasars already existing in our spectroscopic 
quasar-galaxy pair sample, we searched for all neighbouring 
galaxies within the projected distance of 350 kpc and |Az| < 0.03 
from the Photo-z catalogue. We also did the same for our cen- 
tral field galaxies in the spectroscopic field pair sample. As a 
result, we obtained in total 3127 quasar-galaxy pairs and 4214 
field galaxy-galaxy pairs. 

3. Results 

3.1. Morphology classification 

The companion galaxies in the sample underwent a morpholog- 
ical classification in order to make a comparative analysis be- 
tween disk type and spheroidal galaxies. 

The first method attempted was based on a Principal 
Component Analy sis (PCA) of spec troscopic objects using the 
method of Connoll y & Szalavl(ll999h . The spectral classification 
coefficient is generally calculated by using two out of five eigen- 
coefficients that build up the expansion of the eigentemplates 
of a galaxy's spectrum. These eigencoefficients can be extracted 
from the SDSS catalog and are defined as: 

eClass - arctan(-eCoe//2/eCoe//l). (2) 

Early-type galaxies have an eClass val ue -0.35 < eClass < 0, 
while late-types have < eClass < 0.55 (IColdwell & Lambasl 
2006). 

Using the SDSS casjobs interface we obtained the eCoeffi- 
cients eCoeff 1 and eCoefF2 for the 305 quasar companion galax- 
ies. Out of the 305 galaxies 183 were defined by their eClass as 
late-type, 34 as early-type, while the remaining 88 did not fit the 
eClass criteria. 
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Table 1. Visual inspection of 1 83 neighbour galaxies detected by eClass as "late-type". The right column shows the results from the 
visual inspection of these eClass-defined late-type galaxies. The 41 remaining early-type galaxies predicted from eClass, did not 
undergo any visual inspection and are therefore not included in the table. 



eCIass morphology classification 



Galaxy type 


eClass results 


Visual inspection 




(#) 


(#) 


Total 


183 


183 


Late-type 


183 


70 


Early-type or indefinable 





113 
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Table 2. Visual inspection of 172 neighbour galaxies detected by eClass as "late-type" with comparison to Galaxy Zoo results. Out 
of the 275 companion galaxies that had Galaxy Zoo information available, 172 galaxies were predicted to be late- type galaxies by 
the eClass function. The right column shows the results from the visual inspection and Galaxy Zoo-data of these eClass-defined 
late- type galaxies. 



eClass vs Galaxy Zoo morphology classification 



Galaxy type 


eClassresults 


Galaxy Zoo 


Visual inspection 




(#) 


(#) 


(#) 


Total 


172 


172 


172 


Late-type 


172 


30 


25 


Early-type or indefinable 





142 


139 
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A visual morphology inspection by eye of the 183 eClass- 
defined late-type galaxies was also done (see Table 1). Out of 
these only 70 were confirmed as late-type, while the rest ap- 
peared either undefinable or as early-types. 70 late-types con- 
firmed out of 183 eClass-defined galaxies is a success rate of 
38%. Clearly using eClass as a predictor of morphology was a 
failure for this data set. 

A second morphology test using morphologies from the 
Galaxy Zoo project (Lintott et al., 2008, 2010) was compared to 
the eClass results (see Table 2). For the 305 companion galaxies 
275 had been classified in the Galaxy Zoo projecfl From this 
sample of 275 with Galaxy Zoo morphologies, eClass yielded 
172 late-type galaxies and 27 early-type. For the 172 eClass 
defined late-type galaxies 30 were classified by Galaxy Zoo as 
"Spirals", 30 as "Ellipticals" and 112 as "Unidentified." 

Of those 30 galaxies selected by the Galaxy Zoo as late-type 
galaxies, 22 were identified as late-types by visual eye inspec- 
tion, while the remaining 8 were not possible to classify. Of the 
30 galaxies selected by the Galaxy Zoo as early-type galaxies, 
28 were classified by eye as early-type galaxies and two could 
not be classified. This demonstrates that Galaxy Zoo provides a 
relatively good estimate of morphology, especially in compari- 
son with eClass. The clear advantage of the method is likely a 
result of the high internal consistency in the Galaxy Zoo mor- 
phological types. 

3.2. Colors of non-AGN galaxies and quasars 

Rest-frame colors of the stellar populations are easily calculated 
by applying Galactic extinction corrections and K-corrections 
to the five observed SDSS filters (u,g,r,i,z). In our study we 
used u-r color which is supposed to have low noise and is the 
most reliable indicator of the bimodality of the color populations 
(Hiateva et al., 2001). 

The Galactic extinction corrected ("dered") magnitudes were 
taken from the SDSS casjobs database for both quasars and 
galaxies. For the galaxies, the K-correction could also be found 
in the SDSS casjobs database. The K-correction for our quasars 
were calculated assuming a universal power-law spectral energy 
distribution (SED) with the optical flux given by fy = v" and 
mean optical spectral index q'=-0.5. Since most of our quasars 
are at low redshift the K-correction is ve ry small and the spec- 
tral index has changed little in this range ( Kennefick & Bursick'. 
[2008.) . The K-correction for the quasar sample can be calculated 
as 

K{z) = -2.5a logil +z)- 2.5log{l + z) (3) 

The absolute magnitudes of both galaxies and quasars were 
calculated via 

Mat, = m„hs + 5 - SlogiDi) - K(z) (4) 

where D^ is the luminosity distance in parsecs. 

Figure[T]shows how the internal extinction corrected Ue-Re 
coloQ for the companion galaxy changes as a function of the 
projected distance between the quasar and the non-AGN galaxy 
companion. We have chosen to perform binning of the data, and 
calculate the average value of the color in each bin since this 
simplifies the visualization of the results. Little change in color 
as a function of projected distance is observed. 

^ The SDSS casjobs interface allows one to obtain the Galaxy Zoo 
classification for any object in their catalog with a corresponding SDSS 
Object ID. 

* The suffix "e" means "internal extinction corrected" in this context. 



3.3. Surface-densities of non-AGN companions 

The density of the environment around quasars may support a 
merger scenario if small over-densities around the quasars are 
observed. To investigate this the surface densities of galaxies 
around quasars were calculated using 

P = N distance I i{d\f - (d2f)7T (5) 

where N is the number of galaxies in one bin (or "annu- 
lus") with outer radius di and inner radius d2. The companion 
galaxies were divided into 10 bins a 35 kpc each in width. The 
erro r bars were generated v ia bootstrap resampling 100 times 
(see lEfron & TibshiraniL[T98 6) and are plotted in figure |2]with a 
68.3% confidence level, corresponding to 1 cr. 

In general, the relative number of red background galaxies 
at distances further than d > 100 kpc decreases with the smaller 
|Az| used. Perhaps this demonstrates the need to use small |Az| 
cuts when searching for galaxies that have proj ected separ ations 
of d > 100 kpc. Indeed, previous studies (e.g. lYeelll983h have 
shown that galaxies that lie within 100 kpc in projected distance 
have the largest chance to be a physical pair 

Figure |2] shows large overdensites of galaxies around the 
quasars. The overdensity of AGN among the neighbours will be 
described in section 13.51 This has been seen before, but only 
in studies with large contaminations of background and fore- 
ground galaxies (e.g. Serberet al., 2006). Ours is the first study 
with spectroscopic galaxies with well-defined redshifts made on 
such a small-scale environment near quasars (d < 350 kpc) that 
yields similar results. However, the results presented herein are 
much more clear cut. This may give support to quasar forma- 
tion via two scenarios. The first involves the collapse of mate- 
rial into huge dark matter halos where smaller satellite galaxies 
form around the massive quasar simultaneously. The second is 
a merger-driven scenario where massive objects get created via 
hierarchical assembling of smaller parts. In the next we will try 
to separate these two possibilities from each other 

Another curious feature in the surface density plots is the 
minimum of blue galaxies in the bin around 150 kpc. The same 
feature is seen in the AGN surface density (right most column 
of figure|2] A larger sample will be needed to determine whether 
this is a statistical artefact or something physical. 

3.4. Comparison to galaxy-galaxy pairs 

One way of identifying if the surface density enhancements are 
related directly to the nature of the quasars or are present for all 
types of galaxies, is by comparing to the environment of field 
galaxies of same luminosity range . We used a field galaxy sam- 
ple derived from lWav et al.l (1201 ll) as described in Section 1241 

Figure [3] show surface density plots of the field galaxy envi- 
ronment using normal (non-AGN) field galaxies with a similar 
absolute magnitude range to the quasar sample. There is a cer- 
tain skewness towards higher luminosities in the quasar hosts 
if comparing to the field galaxies. We have tried to restrict the 
problems of different luminosity distributions by selecting field 
galaxies within the same luminosity range as the quasars. 

We now would like to know if the increase in the surface 
density of field satellites is as steep as the one of neighbour 
galaxies to quasars. Any increase at short projected separations 
could mean that the overdensity around quasars is steeper than 
around field galaxies and forms in different process. Figure 
shows the ratio of surface densities for satellites of quasars ver- 
sus the field sample. For the sake of the sample sizes we in- 
cluded AGN among the neighbours in this calculation. We see 
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Fig. 1. Uf-Re color of galaxies versus distance to quasar. The left hand side plots contain the Ue - Re colors of the bluer companion 
galaxies where Ue - Re < 2.2. The right hand side plots contain all the non-AGN galaxy companions. Three different redshift 
difference cuts |Az| are pictured for the pairs, top: |Az|=0.001, middle: |Az|=0.006, bottom: |Az|=0.012. 



an increase of the ratio at short projected separations. A poten- 
tial issue here is the difference in the redshift distribution be- 
tween the two samples. The field galaxy-galaxy pairs range be- 
tween 0.03 < z < 0. 1 while the quasar-galaxy pairs range 0.03 < 
z < 0.2. We tried restricting our study by using volume-limited 
subsamples of quasar-neighbour and field-neighbours within this 
redshift range 0.03 < z < 0.1, together with a cut on the lumi- 
nosity magnitude on central field galaxies and quasars on Mr > 
-24. The too small sample size made only the largest cut (|Az|) 
any interesting, confirming a steeper increase in the ratio of blue 
neighbours. 

The fiber collision constraint could eventually lead to an un- 
derestimation of the clustering of blue neighbours near quasars 
at short projected separations. We therefore do a similar plot 
with the two photometric redshift samples of quasar-galaxy and 
field-galaxy pairs where we also select the objects within the 
same absolute magnitude ranges and within the field galaxy red- 
shift range 0.03 < z < 0.1 for a direct comparison. In total, 830 
quasar-galaxy pairs and 2907 field galaxy-galaxy were used in 
this plot. Figure|5]confirms that the clumping of blue neighbours 
might be stronger than observed with the spectroscopic neigh- 
bours. The trend significance in the two different plots (figure 
|4]and figure |5]l is the same (~ 3.5cr), but the photometric red- 
shift samples show a 20-25 % steeper clustering of blue galaxies 
around quasars than does the spectroscopic redshift sample. 

Also the effective r~ 19 limit for spectroscopic targeting bi- 
ases us to miss pairs with a low-luminous neighbour at higher 
redshift. While the magnitude limit would not cause any prob- 
lems in case there would be no correlation between companion 



galaxy luminosity and the distance from quasar, and neither is 
too problematic for the case of an increased amount of massive 
galaxies near the quasar, it could on the other hand create diffi- 
culty in detecting an increased amount of less massive galax- 
ies at short projected separations. A more tolerant magnitude 
limit of spectroscopic targeting could permit us to detect less 
luminous companions. We therefore suspect that the clustering 
of low-luminous, blue galaxies around quasars might be even 
higher than observed at the moment. 

At the end of the day, this revisits an old issue about how the 
occurrence and luminosity of AGN relates to the properties of 
the host galaxies- an issue we thought was clear long ago, but 
apparently is not so clear. 

3.5. AGN fraction in neighbours 

We have a small sample consisting of only 69 quasar- AGN pairs, 
but we are still interested to see if we can get any rough indica- 
tion of how the AGN fraction in those might vary as a function of 
distance. Figure |6]represents how the AGN fraction in the sam- 
ple changes as a function of distance from the quasars. In the 
figure, the AGN fraction in annuli drops off steadily until 100- 
150 kpc, where steeply increases again. Out of the quasar-AGN 
pairs, 14 reside in the closest bin. However, the rough estimation 
might be just random features unless the same can be shown with 
a much larger sample of quasars with neighbours. 

A similar plot for the field-galaxy pairs revealed no changes 
in the AGN fraction as function of distance. One might specu- 
late that there could exist two different processes influencing the 
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Fig. 2. Annular surface densities of different galaxies types around quasars. In the left column: blue galaxies with Ue - Re < 2.2, in 
the second column: red galaxies with Ue-Re > 2.2, in the third column: all galaxies, in the most right column: all AGN neighbours. 
Three different redshift difference cuts |Az| are pictured for the pairs, top: |Az| =0.001, middle: |Az|=0.006, bottom: |Az|=0.012. 
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Fig. 5. Results from the photometric redshift samples. Ratio of 
annular surface densities of galaxies around quasars (pquasars) 
and the annular surface densities of galaxies around central non- 
AGN field galaxies {p field) as function of projected distance. The 
field galaxies were chosen to be inside the same magnitude range 
as the quasars in order to compare objects of same mass range, 
having -24 < Mr <-19. The quasar-galaxy pairs were selected to 
be inside the same redshift range as the field galaxy pairs, 0.03 
<z < 0. 1 . The plot includes AGN among the companion galax- 
ies. In the left column: blue galaxies with Ue-Re < 2.2, in the 
right column: all galaxies. |Az| < 0.03. The ratio is normalized 
for the difference in sample size. 
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Fig. 6. The fraction of AGN companions as a function of dis- 
tance from quasars. Only the largest redshift difference cut 
|Az|=0.012 is used. In total 69 quasar-AGN pairs were detected, 
where 14 reside in the closest bin. 



occurrence of AGN in companion galaxies. While one of these 
processes would be more important at low galaxy densities, the 
other one could be directly associated with mergers or interac- 
tions and therefore only is observed in the closest bin. 
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Fig. 3. Annular surface densities of different types of galaxies around central field galaxies. In the first column: blue galaxies with 
Ue-Re < 2.2, in second column: red galaxies with Ug-Re > 2.2, in third column: all non-AGN galaxies, in the most right column: 
all AGN neighbours. Three different redshift difference cuts |Az| are pictured for the pairs, top: |Az|=0.001, middle: |Az|=0.006, 
bottom: |Az|=0.012. 





3.6. Dust in non-AGN neighbour galaxies 

We are also interested in seeing how the dust content in the 
surrounding galaxies varies as a function of distance from the 
quasar. One standard way of measuring the dust content is by 
using the Ha/H/S emission line flux ratio. Consequently, we can 
only use galaxies with fluxes / > in Ha and HyS. 

Figure Q shows no significant correlation of F(Ha/HyS) with 
distance from the quasar. Unfortunately, not many of our galax- 
ies have H/3 line. 




Distance between objects, kpc 

Fig. 7. Dust extinction: the flux ratio F(lia/iy3 ) in blue galaxies 
with Ue-Re < 2.2. Three different redshift difference cuts |Az| 
are pictured for the pairs, top: |Az|=0.001, middle: |Az|=0.006, 
bottom: |Az|=0.012. The individual errors are unweighted in the 
calculations of mean values and error bars. Only objects with 
emission line fluxes that are greater than three times the mea- 
sured error are included. 



3.7. Star-formation rate (SFR) 

By the mid-1970's several studies (e.g. iBushousel [TQSTh sug- 
gested that tidal interactions and shocks could cause in- 
creased star-formation in galaxies. During mergers large galaxy 
disks could be disrupted by tidal effects, causing large gas 
flows and thus fueling star formation in the central re- 
gions (e.g. Barnes & Hernquistl [l996). More recent simulations 
(iMartig & Bournaudi l2008h have suggested that in mergers be- 
tween two galaxies occuring near a large tidal field (e.g. near 
a group of galaxies) will increase the star-formation rate in the 
merger Recently, many groups have al so observed an incr eased 
star-formation rate during mergers (e.g. lNikolic et al.i 1200 4) and 
also in the neighbourhoods of AGN (e.g. CW06). We calculate 
the star formation rate using the Ha /16563 emission line. In star- 
forming galaxies the Ha emission originates from photoioniza- 
tion by massive, short-lived stars. The advantage of the Ha line 
is that it suffers least from extinction (in the optical) and depends 
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Fig. 4. Results from the spectroscopic redshift samples. Ratio of annular surface densities of galaxies around quasars (pquasars) and 
the annular surface densities of galaxies around central non-AGN field galaxies (p field) as function of projected distance. The field 
galaxies were chosen to be inside the same magnitude range as the quasars in order to compare objects of same mass range, having 
-24 < Mr <-19. The plot includes AGN among the companion galaxies. In the left column: blue galaxies with Ue - Re < 2.2, in 
right column: all galaxies. Three different redshift difference cuts |Az| are pictured for the pairs, top: | Az|=0.001, middle: |Az| =0.006, 
bottom: |Az| =0.012. The ratio is normalized for the difference in sample size. 



least on the metallicity of the ionized gas (e.g. lMoustakas et"an 
I2006h . 

Using the SPS S mea sured Ha emission line the 
iBerpvaU & Ronnbackl (1 19951) calibration was utilized to 
estimate the star formation rate of our galaxy sample: 

L(Ha) ^SFRxl.5l* 10^\W] (6) 

and 

L(Ha) = 4-7tDI y/l^o-h 10"-" (7) 

where SFR is the star formation rate in solar masses per year, 
cr and h are the width and height of the Ha emission line, Dl is 
the luminosity distance in Mpc and the emission line luminosity 
is expressed in Watts. 

As an alternative way of measuring the star formation rate 
we used the [On]- emission line. In SDSS related studies [On] 
/13727 is mainly useful when trying to calculate the star forma- 
tion rate in galaxies at redshifts higher than z=0.4 where the Ha 
emission line is redshifted beyond the red limit of the SDSS 
spectrograph. The [O ii] -emission line is more affected by the 
metallicity and internal dust extinction in galaxies. To reduce 
the reddening and inetallicity effects, we use the K98 calibra- 
tion (iKewlev et alll2004 for calculating the SFR from the [O ii]- 
emission line: 

5^7? = (1.4±0.4)x 10"'^'L([Oii]) (8) 



where the L([0 n]) is given in erg s . 

Figure |8] shows the average star-formation rate in blue com- 
panion galaxies calculated from Ha flux versus distance from 
the quasars (left column), as well as SFR calculated from [O ii] 
(right column). According to the figure |8] SFR via Ha remains 
unaffected as a function of distance, while star formation rate 
from the [On] -line increases in proximity to the quasar. The 
cause of this inconsistency could arise from variations in metal- 
licity or dust extinction as a function of distance from the quasar 
Another possibility may be that a number of AGN galaxies were 
not properly identified in the BPT diagram and therefore create 
a false trend in the [O ii] derived SFR. 

Figure |9] shows the average star-formation rate in all non- 
AGN companion galaxies with detected Ha and [O ii]. The SFRs 
from [O ii] and Ha both decrease as a function of distance from 
the quasar, but the correlation in the SFR from [O ii] is still 
clearer 



3.8. Degree of ionized gas in non-AGN neiglibour galaxies 

The flux ratios F([Oiii]/[On]) and F([Oiii]/Hj0) ai-e useful since 
they provide information on the ionization state of the gas in 
the galaxies. High values could occur both in the presence of an 
AGN, in areas with high star-formation activity and from shock- 
heating of the gas. 
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Fig. 8. Star formation rate in blue companion 
galaxies in solar masses per year. Left col- 
umn shows SFR from from Ha-flux. Right col- 
umn shows SFR calculated from [On] -flux. 
Three different redshift difference cuts are pic- 
tured for the pairs, top: |Az|=0.001, middle: 
|Az|=0.006, bottom: |Az|=0.012. Galaxies with 
SFR>100 and individual errors in SFRj,,or 
>1000 are excluded. The individual errors are 
unweighted in the calculations of mean values 
and error bars. Only objects with emission line 
fluxes that are greater than three times the error 
in the emission line flux are included. 
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Fig. 9. Star formation rate in all non-AGN com- 
panion galaxies in solar masses per year Left 
column shows SFR from from Ha-flux. Right 
column shows SFR calculated from [Oii]-flux. 
Three different redshift difference cuts are pic- 
tured for the pairs, top: |Az|=0.001, middle: 
|Az|=0.006, bottom: |Az|=0.012. Galaxies with 
SFR>100 and individual errors in SFRj„or 
>1000 are excluded. The individual errors are 
unweighted in the calculations of mean values 
and error bars. Only objects with used emission 
line fluxes that are greater than three times the 
error in emission line flux are included. 



The flux ratios are sensitive to metallicity and low metallic- 
ities decreases the cooling efficiency which, in tiirn, in crease s 
the electron temperature and give higher ratios (iSearlel ITqtTI) . 
The ratios may be further modified b y shocks, whose pres- 
ence also will increase the ratios (e.g. iRavmond et al.L Il979t 
ISutherland & DopitaL[T993h . 

In figureflOlthe F([0 iii]/[0 ii])-parameter shows little change 
in ionization for blue galaxy companions as a function of dis- 
tance from the quasar. The same holds true for the F([0 n]fH/3)- 
ratio shown in the same figure. In figure [TT]the number of galax- 
ies with measurable ionization is plotted. The fact that the ion- 
ization cannot be measured for the entire sample makes it dif- 



ficult to drawn any substantive conclusions. We cannot observe 
any correlation of degree of ionization of the gas in the neigh- 
bour galaxies as a function of distance to the quasars. Clearly we 
suffer the fact that only a small fraction of the sample could be 
examined due to the lack of line data. This hinders us to draw 
any firm conclusions. 



4. Discussion and Conclusions 

Using a number of derived quantities from the photometry and 
spectroscopy of the SDSS DR7 we have investigated the proper- 
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Fig. 10. Ionization in blue companion 
galaxies. Left column shows the flux ratio 
F([Oiii]/[Oii]). Right column shows the flux 
ratio F([0 iiil/HjS). Three dilferent redshift 
difference cuts are pictured for the pairs, 
top: |Az|=0.001, middle: |Az|=0.006, bottom: 
|Az|=0.012. Galaxies with ionization>7 are not 
included since they are believed to be AGN 
that have contaminated the sample despite 
using the BPT criteria. Galaxies with individ- 
ual errors > 1000 in ionization are excluded. 
The individual errors are unweighted in the 
calculations of mean values and error. Only 
objects with emission line fluxes greater than 
three times the error in emission line flux are 
utilized. 
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Fig. 11. Left column: number of galaxies with [Oiii] and [On] 
measured. Right column: number of blue galaxies with [O m] 
and H/3 measured. The three |Az| are represented by different 
bar colors. Galaxies with ionization>7 are not included since 
they are believed to be AGN that have contaminated the sam- 
ple despite using the BPT criteria. Galaxies with individual 
errors > 1000 in ionization are excluded. The individual errors are 
unweighted in the calculations of mean values and error Only 
objects with emission line fluxes greater than three times the er- 
ror in emission line flux are utilized. 



ties of 305 quasar-galaxy associations at redshifts of 0.03 < z < 
0.2 and projected distances less than 350 kpc. 

Three different methods for morphology classification of the 
quasar companion galaxies were utilized: eClass, Galaxy Zoo 
and by eye. Galaxy color and spectral line measurements were 
coiTected for internal dust extinction in the non- AGN companion 
galaxies. AGN were separated from th e other neighbour galax- 
ies using BPT-diagram combined with iKauffmann et al.l (l2003h 



criteria and Ha line width. The degree of ionization, dust ex- 
tinction, colors and star-formation rates (SFR) for the non-AGN 
galaxies was measured as function of distance from the quasars. 
We also estimated the surface density distribution of non-AGN 
neighbours and AGN neighbours. Some conclusions are: 

1 . No correlation was found between color, SFR, degree of ion- 
ization or dust extinction in the neighbour galaxy and pro- 
jected distance from the quasar 

2. There appears to be a large increase in the surface density of 
all galaxy types closer than 100 kpc to the quasar. The in- 
crease is larger than for field galaxies of the same luminosity 
range. 

4.1. Star-formation rates and colors 

We see no change in star-formation rate from Ha, color, degree 
of ionization or dust extinction in the companion galaxies as a 
function of projected distance from the quasars. 

This is contrary to what several other studies have reported 
on quasar environments (e.g. CW06) where change in colors and 
SFR in close neighbours have been reported at short projected 
separations. Different studies could have differing SFR conclu- 
sions because the methods used to separate AGN from normal 
galaxies vary between the studies. 

The general lack of correlations in the quasar non-AGN com- 
panions as function of distance from the quasar could be a result 
from having a sample consisting mainly of faint quasars. If one 
wishes to see effects from the hypothesized AGN feedback on 
the distance scale up to 350 kpc, one might need to perform 
a study with a sample consisting of very strong high-redshift 
quasars. Since the present day catalogues of quasars have too 
few observable companion galaxies we cannot do a similar high- 
redshift study yet. At the moment the SDSS DR7 does not allow 
either study due to the far too small sample size one obtains from 
the present sample and the lack of objects with spectroscopic 
redshifts above z > 0.2. 
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The lack of observable effects on companions can also be 
a consequence of quasar influence being episodic rather than 
continuous. A short episode of effects from the quasar on the 
companion could perhaps not be disentangled by the color, star- 
formation and ionization tests unless the gas would be com- 
pletely swept away. Ionization from the quasar on the interstellar 
medium in the companion galaxy might regenerate to a normal 
state with time and a past burst of star-formation in the near- 
est neighbours could already have died out, leaving no trace in 
our observations. Therefore it might be more advantageous for 
future studies to concentrate on a combination of total masses 
of the atomic and molecular gas content together with spectral 
history modeling for the nearest neighbours to disentangle the 
possibilities of a past episode of starburst or quenching. 

We could still expect to see effects of clos e interactions on 
the SFR or color as many other studie s (e.g. iKennicutt et al.L 
Il987t lEmsonetalL 120101: iNikolic et al.L 12004 reported obser- 
vations that interactions and mergers can enhance the SFR. 
However there are contraindications that interactions are not 
enough for causing starbursts unless the galaxies actually are 
merging, and that even mergers might not be a sufficient con- 
dition (e.g. lBergvall et al.l l2003l) . This might be consistent with 
that we cannot see any increase of SFR or change in color in our 
quasar neighbours at short projected separations. 

4.2. AGN fraction of neighbours 

Our results suggest that there might be a drop of the annular 
AGN fraction among the neighbours to quasars, and a steep in- 
crease the closest 50-100 kpc. 

A larger study of AGN pairs by Liu et al. (2011) shows 
steadily dropping cumulative AGN fraction represents a large 
sample with 1286 AGN- AGN pairs. Their pairs are a mixture of 
all kind of active galaxies with a S/N ratio > 3 for each emission- 
line used in the AGN identification; narrow- line galaxies derive d 
from the MPA-JHU S PSS DR7 ca t alogue (lAiharaet al.Ll201lh . 
narrow-Une quasars jReves et all l2008l) . broad-line galaxies 
from (ffaoetaT, '2005") and broad-line quasars from the quasar 
catalogue (Schneider et al., 2010). In Liu's study they show that 
for AGN-pairs with emission lines having a S/N ratio>15, 80% 
of the pairs display tidal features, something that supports merg- 
ing and tidal interactions behind AGN activity. 

The implicit assumption for the selection of Liu's parent 
AGN sample is that all AGN are the same type of object ob- 
served with different viewing angle. If there are clearly different 
formation processes for different types of active galaxies and the 
quasars in our sample are not formed in the same way as the 
AGN in their sample, we would see a different behavior of the 
AGN fraction versus projected distance between the objects in 
comparison to our study with quasar-AGN pairs. 

We morphologically examined the 14 quasar-AGN pairs 
within 50 kpc, where 8 showed tidal features. An interesting 
question is whether our they are real AGN or possibly normal 
galaxies contaminated by light from the quasar It is possible 
that if using a higher S/N-cut could decrease the fraction in the 
nearest bin, where the contamination from the quasars also has 
the highest probability to influence our emission lines. Maybe 
it could also be explained in the light echo from quasar such 
that h as been demonstrated in Hanny's Voorwerp (Lintott et al.L 
l2009h . In such case the contamination would be separated 
largely by the time of light travel between the two quasar and 
the neighbour, in such a way that the quasar already would 
have evolved. We note that these 14 neighbours have rather low 
F([0 ra]/Hy6) ~ 1 . X-ray measurements indicative of non-thermal 



source could be a good way to find out if they are real AGN or 
not. 



4.3. Enfianced surface densities 

We plotted the surface densities of the blue neighbour galaxies, 
the red galaxies and all galaxies combined and found a large in- 
crease in surface density of all galaxy types at bins closer than 
150 kpc to the quasar. Especially high is the increase in blue 
gas-rich galaxies and AGN if one compares to the background 
of galaxies further away from the quasar (d > 150 kpc). In case 
mergers are responsible for quasar activity, one could imagine 
that the nearby blue galaxies and AGN also are products of local 
mergers. But this hypothesis would require modeling of star for- 
mation history of the objects, which is beyond the scope of this 
paper. 

At the same time an increased density around the quasars 
does not necessarily imply a merger scenario, even though the 
temptation to assign this explanation to the result might be 
strong. Monolithic collapse models can predict the formation of 
smaller objects around more massive ones when both smaller 
and larger objects collapse in the gr avitational well of the same 
halo (e.g. iPress & Schechteil 11974 . The ACDM model itself 
predicts a large number of satellite galaxies around massive 
galaxies, something that is still difficult to know whether it is 
the case or not here. 

Our comparison to field galaxies of the same luminosity dis- 
tribution as the quasars shows that at short projected separations 
quasars might have more satellites than the field galaxies of sim- 
ilar luminosity distribution, and particularly blue satellites. This 
is possible to interpret as that quasars form from mergers of gas- 
rich (blue) objects and do not form in the same way as field 
galaxies. 

However, the most interesting feature of the data might lie in 
the difference between the total number of non-AGN neighbours 
to quasars or field galaxies in the |Az|=0.001 and |Az|=0.012 
samples, see table 3. As indicated from the table, the fraction 
of satellites residing in the smaller subsample is larger around 
quasars than around field galaxies, which is another support for 
stronger clustering of satellites around the quasars. 

Also CW06 found support in their studies of the surface 
density distributions of galaxies close to clusters, close to field 
galaxies and close to quasars that galaxies in the neighbourhood 
of field galaxies are found in lower density environment than 
galaxies around q uasars. Similar results have been obtained by 
other groups (e.g. Serber et al.i l2006i IStrand et al.l I2008D . 

How could we further test whether the increased number of 
satellites around quasars is a product of hierarchical structure 
formation or of monolithic collapse? In a merger scenario, the 
observed neighbour galaxies were not involved in forming our 
quasars, and we see only those satellites that are left after the 
process. One way would be to construct a SDSS sample of field 
quasars at low-redshifts (0.03 < z < 0.2) but that have no com- 
panion closer than 350 kpc. In case of merger scenario being the 
secret of quasar activity, these quasars should have an average 
luminosity that is fainter than in the sample with quasars that 
have companion galaxies. A second way would be to use the lo- 
cal velocity dispersions of neighbours as indicator of tendency 
to merging. However, this would demand much larger samples 
as the lack of accuracy in SDSS redshifts might make it difficult 
to draw any conclusions otherwise. 

An additional way is to approach the problem is by using 
modelling to predict the number of lower-mass objects around 
quasars of the average luminosity in our quasar sample (e.g. 



13 



Beatriz Villarroel: In the Neighbourhood of Tame Monsters. 



Table 3. Number of non-AGN satellites around quasars and field galaxies within 350 kpc. The number of satellites in the subsample 
|Az|=0.001 (column 2) is compared to the main sample |Az| =0.012 (column 3). The fraction of satellites in the subsample |Az|=0.001 
(column 4). 



Clustering of neighbours 



Central galaxy 


|Az|=0.001 


|Az|=0.012 


fraction 




(#) 


(#) 


(dimensionless) 


Quasars 


138 


235 


59% 


Field 


161 


326 


49 % 



ICarlber4ll99d ). Similar attempts of investigating the clustering 
of quasars have been done (e.g. Shen et al., 2010) where an ex- 
cess of clustering could be shown, but where the results could 
not be clearly put in context of any given structure formation 
scenario. 

4.4. Missing AGN and blue galaxies? 

In this section we would like to offer the reader some more spec- 
ulative statements on the data aiming to nurture a curious mind, 
yet without the intention of being conclusive or complete. 

A curious thing about the surface density plots is the mini- 
mum of blue galaxies and AGN in the bins around ~ 150 kpc 
on a statistical significance level of ~ 1.5 cr. This could be a 
product of noisy data or low number statistics, but at the same 
time we see a slight increase of red galaxies in the 35 kpc closer 
surface density bin. We did similar surface density plots for the 
photometric pair samples with only color separation (blue and 
red) among the neighbours, see figure [12] For the quasar-galaxy 
pairs with a blue neighbour (1732 pairs), the 150 kpc gap again 
could be spotted on a ~ 2 cr significance level. Nothing simi- 
lar was spotted for the photometric field galaxy pairs with blue 
companions. 

Could this mean that those missing blue galaxies got trapped 
by the gravitational field of the quasar and got their gas stripped 
off (hence became red) as they fell into circular orbits around the 
quasar or lost it due to ionization pressure from the quasars ? 

It should be noted that the surface densities are shown on a 
logarithmic scale. It may seem like the loss of blue galaxies and 
AGN neighbours in the ~ 150 kpc bin is larger than the gain 
of red galaxies. Perhaps some galaxies lost the gas when falling 
into bound orbits. Or maybe there is an unseen population of 
very faint galaxies in this bin. 

What if the galaxies simply not had a chance even to form 
any stars in their recent star formation history, before get- 
ting trapped b y the quasars? A study of Milky Way satellites 
dGrcevich & P utman. 2009) has revealed that within 270 kpc of 
the Milky Way, almost all satellite galaxies have lost their gas 
and are spheroidal. This is support for tidal stripping processes 
within the distance range in our own study. What if something 
similar has happened to our the blue non-AGN companions at 
~ 150 kpc? These highly speculative scenarios would however 
need careful measurements of the velocity fields of the close 
quasar companion galaxies in order to follow how the dynam- 
ics of the gaseous content of the companion galaxies changes 
around the minimum. 

Many thoughts remain as pure speculations until we find a 
way to increase our neighbour galaxy sample without losing 
precision in redshift estimates. It is difficult to know whether 
the minimum is real or simply an artefact from having too few 
quasar-galaxy pairs. Further studies with and increased number 
of galaxies will be required to confirm (or reject) this detection. 
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Fig. 12. Results from the photometric redshift 
samples. Annular surface densities of neigh- 
bour galaxies around quasars. In the first col- 
umn: all blue neighbours with Ue - Re < 2.2, 
in middle column: all red neighbours with (/<, - 
Re > 2.2, in the most right column: all neigh- 
bours. The total number of galaxies is indicated 
in each plot. Only one redshift difference cut 
|Az| < 0.03 is used. AGN are included among 
the neighbours. 
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